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We present the first measurement of tfie inclusive three-jet differential cross section as a function 
of the invariant mass of the three jets with the largest transverse momenta in an event in pp 
collisions at y/s — 1.96 TeV. The measurement is made in different rapidity regions and for different 
jet transverse momentum requirements and is based on a data set corresponding to an integrated 
luminosity of 0.7 fb~^ collected with the DO detector at the Fermilab Tevatron Collider. The results 
are used to test the three-jet matrix elements in perturbative QCD calculations at next-to-leading 
order in the strong coupling constant. The data allow discrimination between parametrizations of 
the parton distribution functions of the proton. 

PACS numbers: 13.87.Ce, 12.38.Qk 



The production cross section for jets with large trans- 
verse momenta (px) with respect to the beam axis in 
hadron-hadron collisions is predicted by perturbative 
QCD (pQCD) and is sensitive to the strong coupling con- 
stant (as) and the parton distribution functions (PDFs) 
of the hadrons. Deviations from the pQCD predictions 
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may indicate the presence of physics processes not in- 
cluded in the standard model. Recent measurements of 
inclusive jet and dijet production in pp collisions at a 
center-of-mass energy of ^/s — 1.96 TeV [ll-Q have been 
used to determine as [7\ and the proton PDFs and 
to set limits on a number of models of physics beyond the 
standard model [53 . This demonstrates the success of 
pQCD in describing observables which are directly sensi- 
tive to the matrix elements of 0{al). Testing pQCD at 
higher orders of requires measuring cross sections for 
higher jet multiplicities. 

The three-jet cross section is directly sensitive to the 
pQCD matrix elements of 0{a'^), and therefore has a 
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FIG. 1: (Color online.) The differential cross section dazjct / dMy^^t (a) in different rapidity regions and (b) for different 
PT3 requirements. The solid lines represent the NLO pQCD matrix element calculations using MSTW2008NLO PDFs and 
as{Mz) ~ 0.1202 which are corrected for non-perturbative effects. 



higher sensitivity to as compared to inclusive jet and 
dijet cross sections, while having a similar sensitivity to 
the PDFs. Since pQCD calculations are available to next- 
to- leading order (NLO) in as lll-ll4l|. the three-jet cross 



section can be used for precision phenomenology such as 
simultaneous determinations of as and PDFs from ex- 
perimental data. In such QCD analyses [1, 15|, the in- 
formation from three-jet cross sections can supplement 
that from inclusive jet and dijet cross sections, partially 
decorrelating the results for and the PDFs. 

In this Letter, we present the first measurement 
of the inclusive three-jet differential cross section, 
dcTsjet/dMsjet, in Pp collisions at y/s — 1.96 TeV, as 
a function of the invariant mass (Majet) of the three 
highest-pT jets in each event. The data sample, collected 
with the DO detector ^ during 2004-2005 in Run II of 
the Fermilab Tevatron Collider, corresponds to an inte- 
grated luminosity of 0.7 fb~^. In the experiment and in 
the theoretical calculations used in this analysis, jets are 
defined by the Run II midpoint cone jet algorithm 17 1 
with a cone of radius 7?,conc = 0.7 in rapidity y and az- 
imuthal angle cj). Rapidity is related to the polar scat- 
tering angle 6 with respect to the proton beam axis by 
y = i In [(1 -I- /3 cos0)/(l — /3cos0)], where /? is defined 
as the ratio between momentum and energy (/3 = \p\/E). 
The inclusive three-jet event sample consists of all events 
with three or more jets which pass given px and \y\ re- 
quirements. The Msjet dependence of the inclusive three- 
jet cross section is measured for five scenarios with dif- 
ferent jet pt requirements and in different regions of jet 
rapidity. Jets are ordered in descending px and the pr re- 
quirements are pri > 150 GeV and prs > 40 GeV (with 
no further requirement for pt2)- The rapidities of the 



three leading pr jets are restricted to \y\ < 0.8, \y\ < 1.6, 
or \y\ < 2.4, in three different measurements. Two ad- 
ditional measurements are made for pTs > 70 GeV and 
Pt3 > 100 GeV, both requiring \y\ < 2.4. For jets defined 
by the cone radius 7?.conc and for a given pxa requirement, 
the relative pt between two jets {kj_) could be as low as 
kj_ ~ T^conc ■ PT3, which introduces an additional, softer 
scale in the process (since k± < pT3 for 7?.conc — 0.7). 
The phase space with k± below the pxa requirement can 
be avoided by an additional requirement on the angular 
separation of the three leading px jets. In all scenarios, 
all pairs of the three leading px jets are required to be 
separated by AR = ^/{Ay)^ + (A0)2 > 1.4 (= 2•7^cone)■ 
With this separation requirement, the smallest accessi- 
ble fc_L of the jets is always above prs- Furthermore, 
this separation requirement also reduces the phase space 
in which pairs of the three leading px jets are subject 
to the overlap treatment in the cone jet algorithm jl7| . 
Since the overlap treatment can strongly depend on de- 
tails of the energy distributions in the overlap area, this 
region of phase space may not be well modeled by pQCD 
calculations at lower orders. In the remaining analysis 
phase space, NLO pQCD calculations are not affected by 
the Run II cone algorithm's infrared sensitivity [18]. The 
data are corrected for instrumental effects and are pre- 
sented at the "particle level," which includes all stable 
particles as defined in Ref. |19| . 

A detailed description of the DO detector can be found 
in Ref. [l^. The event selection, jet reconstruction, and 
jet energy and momentum correction in this measure- 
ment follow closely those used in our recent inclusive jet 
and dijet measurements Jets are reconstructed in 

the finely segmented DO liquid-argon/uranium calorime- 
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FIG. 2: (Color online.) Ratios of the differential cross sections dcrsjet /dMsjet measured in different rapidity regions and for 
different pra requirements and the pQCD predictions for different PDFs. The inner uncertainty bars indicate the statistical 
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ter which covers most of the solid angle for polar angles 
of 1.7° S < 178.3° For this measurement, events 
are triggered by the jet with highest pr- Trigger efficien- 
cies are studied by comparing the three-jet cross section 
in data sets obtained by inclusive jet triggers with differ- 
ent pt thresholds in regions where the trigger with lower 
threshold is fully efficient. The trigger with lowest pt 
threshold is shown to be fully efficient by studying an 
event sample obtained independently with a muon trig- 
ger. In each Majet bin, events are taken from a single 
trigger which is chosen such that its efficiency is above 
99%. 

The position of the pp interaction is determined from 
the tracks reconstructed in the silicon detector and scin- 
tillating fiber tracker located inside a 2 T solenoidal mag- 
net [16i] . The position is required to be within 50 cm 
of the detector center in the coordinate along the beam 
axis, with at least three tracks pointing to it. These 
requirements discard (7.1-8.6)% of the events. Contri- 
butions from cosmic ray events are suppressed by requir- 
ing the missing transverse energy (-^t) in an event to be 
< 0.5 ■pTi- This requirement is applied before the jet 
four-momenta are corrected, and its efficiency for signal 
is found to be > 99.5% Requirements on character- 
istics of calorimeter shower shapes are used to suppress 
the remaining background due to electrons, photons, and 
detector noise that would otherwise mimic jets. The effi- 
ciency for the shower shape requirements is above 97.5%, 
and the fraction of background events is below 0.1% for 



all M: 



The jet four-momenta reconstructed from calorimeter 
energy depositions are then corrected, on average, for the 
response of the calorimeter, the net energy flow through 
the jet cone, additional energy from previous and consec- 
utive beam crossings, and multiple pp interactions in the 
same event, but not for muons and neutrinos. The abso- 
lute energy calibration is determined from Z — > ee events 
and the pt imbalance in 7 -I- jet events in the region 
\y\ < 0.4. The extension to larger rapidities is derived 
from dijet events using a similar data-driven method. In 
addition, corrections are applied which take into account 
the difference in calorimeter response due to the differ- 
ence in the fractional contributions of quark- and gluon- 
initiated jets in the dijet and the 7 -t- jet event samples. 
These corrections of the order (2-4)% are determined us- 
ing simulated jets produced with the PYTHIA event gener- 
ator [2^ that have been passed through a GEANT-based 
detector simulation [2l|. The total corrections for the 



jet four-momenta vary between 50% and 20% for jet pr 
between 50 and 400 GeV. These corrections adjust the 
reconstructed jet energy to the energy corresponding to 
the stable particles that entered the calorimeter except 
for muons and neutrinos, which are accounted for later 
by a separate correction. An additional correction is ap- 
plied for systematic shifts in\y \ due to detector effects [J]. 
The three-jet invariant mass is then computed from the 
corrected jet four-momenta of the three highest-pT jets 
in the event. 
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FIG. 3: (Color online.) Ratios of the differential cross sections dcrsjet /dMsjet measured in different rapidity regions and for 
different pra requirements and the pQCD predictions for different PDFs. The inner uncertainty bars indicate the statistical 
uncertainties, and the total uncertainty bars display the quadratic sum of the statistical and systematic uncertainties. 



The differential cross sections da^\et I dM^Sat are cor- 



rected for experimental effects [22|. Particle-level 
jets from events generated with SHERPA [ISl with 
MSTW20Q8LO PDFs % are processed by a fast simu- 
lation of the DO detector response. The simulation is 
based on parametrizations of resolution effects in p^, the 
polar and azimuthal angles of jets, jet reconstruction effi- 
ciencies, and misidentification of the event vertex, which 
are determined either from data or from a detailed simu- 
lation of the DO detector using GEANT. The pt resolution 
for jets is about 15% at 40 GeV, decreasing to less than 
10% at 400 GeV. The generated events are reweighted to 
match the Majct, pr, and |y| distributions in data. To 
minimize migrations between Majot bins due to resolu- 
tion effects, we use the simulation to obtain a rescaling 
function in reconstructed Afajct that optimizes the cor- 
relation between the reconstructed and true values. The 
bin sizes in the Majot distributions are chosen to be ap- 
proximately twice the Afajet resolution. The bin purity 
after A^ajct rescaling, defined as the fraction of all recon- 
structed events that were generated in the same bin, is 
above 40% for all bins. We then use the simulation to 
determine Majot bin correction factors for instrumental 
effects for the differential cross sections in the five differ- 
ent scenarios. These also include corrections for the en- 
ergies of unreconstructed muons and neutrinos inside the 
jets. The total correction factors for the differential cross 
sections vary from about 1.0 at Afajct = 0.4 TeV to 1.1 at 
1.1 TeV for \y\ < 0.8 and between 0.89 at Mgjot = 0.4 TeV 
to 0.96 at 1.1 TeV for \y\ < 2.4. The dependence of the 
correction factors on the reweighting function is taken 



into account as an uncertainty. The corrected differential 
cross section in each scenario is presented at the "particle 
level" as defined in Ref. [IS]. 

In total, 65 independent sources of experimental sys- 
tematic uncertainties are identified, mostly related to jet 
energy and jet pr resolution. The effects of each source 
are taken as fully correlated between all data points. The 
dominant uncertainties for the differential cross sections 
are due to the jet energy calibration [(10-30)%] , the lumi- 
nosity uncertainty (6.1%), and the jet pT resolution [(1- 
5)%]. Smaller contributions come from the uncertainties 
in systematic shifts in y (3%), reweighting of the gener- 
ated events (2.5%), trigger efficiency uncertainties (2%), 
and from the jet 9 resolution (1%). All other sources are 
negligible. The systematic uncertainties are never larger 
than 30%, and for Majet < 0.9 TeV, they are between 
11% and 20%. 

The results for the differential cross sections for differ- 
ent rapidity and pt3 requirements are given in Table |T] 
and displayed in Fig. [T] A detailed documentation of the 
results, including the individual contributions from all 
65 sources of correlated uncertainties is provided in the 
supplemental material . The quoted central values of 
Msjet at which the data points are presented are the loca- 
tions where the bin averages have the same value as the 
differential cross section [25[ , as determined using smooth 
parametrizations of the data. The data are compared to 
theory predictions which have been obtained from NLO 
pQCD calculations with non-perturbative corrections ap- 
plied. The non-perturbative corrections are determined 
using PYTHIA with "tune DW" ^2d\. They are defined 
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TABLE I: The thrcc-jct differential cross section dasjet /rfAfsjet and the theoretical predictions based on NLO pQCD (for 
MSTW200NLO PDFs with as{Mz) = 0.1202) plus non-perturbative corrections, for renormalization and factorization scales 

l^r = Hf = (PTI +PT2 +PT3)/3. 



Mass 


Central 


Measured 


Statistical Systematic 


Theory 


Non-perturbative corrections 




range 


mass 


cross section 


uncertainty uncertainty 


cross section 


Hadronization 


Underlying 


Total 


TeV 


TcV 


pb/TeV 


% 


% 


pb/TcV 


% 


event % 


% 


for \y\ < 0.8 and prs > 40 GeV 


0.40-0.45 


0.424 


1.93 X 


lO'' 


1.8 


+11.2, -11.5 


1.99 X lO'' 


-9.0 


9.1 


-0.7 


0.45-0.50 


0.474 


1.12 X 


10^ 


2.2 


+12.0, -11.1 


1.17 X 10^ 


-8.0 


9.0 


-0.2 


0.50-0.55 


0.524 


5.89 X 


10^ 


2.9 


+ 12.2, -12.1 


6.08 X lOi 


-7.5 


8.7 


0.6 


0.55 0.61 


0.578 


2.76 X 


10^ 


3.8 


+ 11.6, -12.0 


2.97 X lOi 


-7.1 


8.2 


0.5 


0.61 0.67 


0.638 


1.22 X 


10^ 


5.6 


+ 12.5, -11.8 


1.33 X lOi 


-6.9 


8.0 


0.6 


0.67-0.74 


0.703 


5.75 X 


10° 


6.4 


+ 12.4, -13.5 


5.98 X 10° 


-6.7 


8.0 


0.7 


0.74-0.81 


0.773 


2.56 X 


10° 


9.6 


+14.5, -13.6 


2.45 X 10° 


-6.8 


7.5 


0.2 


0.81-0.90 


0.851 


8.47 X 


10-1 


14.8 


+14.9, -14.3 


8.99 X 10-1 


-6.9 


6.3 


-1.1 


0.90-1.10 


0.976 


1.85 X 


10-1 


21.3 


+18.1, -17.2 


1.64 X 10-1 


-7.0 


5.3 


-2.1 


for \y\ < 1.6 and prs > 40 GeV 


0.40-0.45 


0.434 


1.01 X 


10^ 


0.8 


+11.7, -12.3 


1.10 X 10^ 


-10.0 


8.4 


-2.4 


0.45-0.50 


0.476 


8.74 X 


10^ 


0.9 


+13.3, -12.5 


9.95 X 10^ 


-8.9 


9.9 


0.2 


0.50-0.55 


0.525 


6.28 X 


10^ 


1.0 


+12.6, -13.0 


7.18 X 10^ 


-8.3 


10.6 


1.4 


0.55-0.61 


0.579 


3.95 X 


10^ 


1.1 


+14.3, -12.5 


4.47 X 10^ 


-8.2 


10.7 


1.7 


0.61 0.67 


0.639 


2.08 X 


10^ 


1.5 


+14.5, -13.9 


2.42 X 10^ 


-8.0 


10.5 


1.6 


0.67-0.74 


0.703 


1.00 X 


10^ 


1.9 


+14.7, -14.9 


1.18 X 10^ 


-7.8 


10.1 


1.5 


0.74-0.81 


0.773 


4.54 X 


lOi 


2.8 


+16.4, -15.9 


5.14 X lOi 


-7.7 


9.8 


1.4 


0.81-0.90 


0.851 


1.64 X 


10^ 


4.0 


+17.2, -15.6 


1.91 X lOi 


-7.8 


9.8 


1.2 


0.90-1.10 


0.978 


3.18 X 


10° 


6.0 


+20.3, -20.8 


3.47 X 10° 


-8.2 


9.5 


0.5 


1.10-1.50 


1.215 


8.71 X 


10-2 


20.6 


+27.4, -26.5 


1.03 X 10-1 


-9.1 


8.4 


-1.5 


for \y\ < 2.4 and pra > 40 GeV 


0.40-0.45 


0.419 


1.23 X 


lo-" 


0.7 


+ 12.1, -12.5 


1.33 X lO'' 


-10.5 


8.6 


-2.8 


0.45 0.50 


0.477 


1.17 X 


10=* 


0.8 


+13.6, -12.5 


1.30 X 10^ 


-9.6 


9.9 


-0.7 


0.50-0.55 


0.526 


9.23 X 


10^ 


0.8 


+13.0, -13.1 


1.05 X 10^ 


-9.2 


10.5 


0.3 


0.55-0.61 


0.580 


6.68 X 


10^ 


0.9 


+14.6, -13.3 


7.47 X 10^ 


-9.1 


10.7 


0.7 


0.61-0.67 


0.639 


4.21 X 


10^ 


1.1 


+14.9, -14.7 


4.81 X 10^ 


-9.1 


10.8 


0.8 


0.67-0.74 


0.704 


2.48 X 


10^ 


1.3 


+15.5, -15.8 


2.81 X 10^ 


-9.2 


11.0 


0.8 


0.74-0.81 


0.773 


1.33 X 


10^ 


1.7 


+17.3, -15.6 


1.50 X 10^ 


-9.4 


11.6 


1.0 


0.81-0.90 


0.852 


6.01 X 


10^ 


2.3 


+18.1, -18.9 


6.96 X lOi 


-9.8 


12.3 


1.3 


0.90-1.10 


0.983 


1.38 X 


10^ 


3.1 


+20.6, -21.2 


1.60 X IQi 


-10.0 


12.8 


1.5 


1.10-1.50 


1.225 


5.40 X 


10-1 


10.3 


+29.5, -27.2 


6.54 X IQ-i 


-10.1 


13.7 


2.2 


for \y\ < 2.4 and prs > 70 GeV 


0.40-0.45 


0.421 


3.84 X 


10^ 


1.4 


+ 12.2, -12.5 


4.24 X 10^ 


-11.6 


7.5 


-5.0 


0.45-0.50 


0.481 


4.29 X 


10^ 


1.3 


+13.0, -12.7 


4.91 X 10^ 


-10.1 


9.2 


-1.8 


0.50-0.55 


0.526 


3.71 X 


10^ 


1.3 


+13.2, -13.6 


4.29 X 10^ 


-9.5 


10.0 


-0.5 


0.55-0.61 


0.580 


2.87 X 


10^ 


1.3 


+13.7, -13.5 


3.29 X 10^ 


-9.5 


10.1 


-0.3 


0.61-0.67 


0.639 


1.90 X 


10^ 


1.6 


+14.8, -14.0 


2.22 X 10^ 


-9.5 


10.1 


-0.3 


0.67-0.74 


0.704 


1.20 X 


10^ 


1.9 


+15.6, -15.5 


1.38 X 10^ 


-9.5 


10.3 


-0.2 


0.74-0.81 


0.774 


6.52 X 


lOi 


2.5 


+16.7, -14.8 


7.68 X IQi 


-9.6 


10.8 


0.2 


0.81-0.90 


0.853 


3.18 X 


lOi 


3.1 


+17.6, -19.6 


3.80 X IQi 


-9.9 


11.6 


0.6 


0.90-1.10 


0.985 


8.27 X 


10° 


4.0 


+21.5, -20.7 


9.71 X 10° 


-10.2 


12.1 


0.7 


1.10-1.50 


1.235 


3.40 X 


10-1 


13.0 


+28.6, -24.5 


4.21 X 10-1 


-10.3 


12.7 


1.1 












for \y\ < 2.4 and pra '. 


> 100 GcV 








0.40-0.45 


0.424 


3.80 X 


lOi 


4.0 


+13.9, -13.5 


4.80 X lOi 


-13.8 


5.3 


-9.3 


0.45-0.50 


0.472 


7.55 X 


lOi 


2.8 


+15.0, -11.9 


8.72 X lOi 


-11.2 


6.6 


-5.3 


0.50-0.55 


0.535 


8.53 X 


lOi 


2.6 


+13.6, -13.7 


9.61 X lOi 


-9.9 


7.7 


-2.9 


0.55-0.61 


0.581 


7.56 X 


lOi 


2.4 


+12.9, -13.7 


8.66 X lOi 


-9.5 


8.6 


-1.7 


0.61-0.67 


0.640 


5.59 X 


lOi 


2.8 


+13.9, -13.1 


6.60 X lOi 


-9.2 


9.3 


-0.7 


0.67-0.74 


0.704 


3.84 X 


lOi 


3.1 


+15.0, -14.7 


4.46 X lOi 


-8.9 


9.8 


0.1 


0.74-0.81 


0.774 


2.18 X 


lOi 


4.0 


+18.9, -16.4 


2.69 X lOi 


-8.8 


10.0 


0.3 


0.81-0.90 


0.853 


1.16 X 


lOi 


4.8 


+16.5, -18.6 


1.41 X lOi 


-9.3 


10.1 


-0.1 


0.90-1.10 


0.985 


3.38 X 


10° 


6.0 


+22.2, -20.7 


3.89 X 10° 


-9.8 


10.4 


-0.5 


1.10-1.50 


1.225 


1.58 X 


10-1 


18.4 


+25.0, -22.2 


2.10 X 10-1 


-10.1 


11.5 


0.3 
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as the combination of the corrections due to hadroniza- 
tion and underlying event and vary between —10% and 
+2% (given in Table Using different pythia set- 
tings (A, BW, Zl, Perugia soft, Perugia hard tunes) af- 
fects the individual corrections by less than half of their 
sizes and the total corrections by less than 5%. The 
NLO results are c omp uted using fastnlo ^ based 
on NLOJET-h-h O [3 with MSTW2008NLO PDFs Q 
and the corresponding value of as{Mz) — 0.1202. The 
central choice /xg for the renormalization and factoriza- 
tion scales is the average pr of the three leading pT jets 
= ^/ = ^0 = {PTi + PT2 + Pts)/^- For a direct 
comparison of the theoretical predictions with data, the 
ratio of data and theory is displayed by the markers in 
Fig. [2] for all five scenarios. The effects of independent 
variations of renormalization and factorization scales be- 
tween /xo/2 and 2/io are displayed by the dotted lines. 
These variations affect the predicted cross sections be- 
tween + {5-10)% and -(15-20)%. 

The MSTW2008NLO PDF uncertainties (correspond- 
ing to the 68% C.L.) are shown by the light band. The 
ratios of data and theory are almost constant, with only 
a small dependence on Msjet and the \y\ and prs require- 
ments. The central data values are below the central 
theory predictions, by approximately (4-15)% in the dif- 
ferent scenarios, slightly increasing with \y\ and with pTs- 
In all cases, the data lie inside the range covered by the 
scale variation. 

In addition to the MSTW2008NLO PDFs, Fig.[2]shows 
also predictions for CTIO PDFs and the correspond- 
ing value of as{Mz) = 0.118, normalized by the pre- 
dictions for MSTW2008NLO and represented by the 
solid lines. To compare the CTIO PDF uncertainties 
(which have been published at the 90% C.L.) with the 
experimental uncertainties (corresponding to one stan- 
dard deviation), the former have been scaled by a fac- 
tor of 1/1.645 [ii]. The resulting 68% C.L. uncertain- 
ties are displayed around the CTIO central values by 
the dark band. The CTIO PDFs predict a different 
shape for the M^jct dependence of the cross section. For 
Majot < 0.6 TeV, the central resuhs for CTIO PDFs 
agree with those for MSTW2008NLO, while the CTIO 
predictions at M^jct = 1.2 TeV are up to 30% higher. 
These discrepancies at highest Majet are larger than the 
combined 68% C.L. uncertainty bands of the CTIO and 
MSTW2008NLO PDFs. 

Calculations for additional PDFs are compared to the 
data in Fig. [31 These are the PDF parametrizations 
NNPDFv2.1 10] {ot,{Mz) = 0.119), ABKM09NLO [11] 
(as{Mz) = 0.1179), and HERAPDFvl.O t29|] {as[Mz) = 
0.1176). The results for NNPDFv2.1 agree everywhere 
within ±4% with those from MSTW2008NLO. The cross 
sections predicted for HERAPDFvl.O are (15-20)% be- 
low those for CTIO everywhere and their Afajct distri- 
butions have a similar shape. The A/ajct dependence 
of the calculations for the ABKM09NLO PDFs is be- 



tween the shapes of MSTW2008NLO/NNPDFv2.1 and 
CTlO/HERAPDFvl.O. At low Mgjet, the predictions 
for ABKM09NLO agree with those for HERAPDFvl.O, 
while at higher Majot, they predict the smallest cross sec- 
tions of all PDFs under study. 

The level of agreement between theory and data can 
not be directly judged from the comparisons in Figs. [2] 
and[3l but requires taking into account correlations of ex- 
perimental uncertainties. While some experimental un- 
certainties (like the luminosity uncertainty) allow to shift 
the data points coherently up or down, other uncertainty 
sources (such as the jet energy calibration), have Msjct- 
dependent effects which also allow changes to shapes of 
the data distributions. To quantify the significance of 
the differences between theory and data in normalization 
and shape as observed in Figs. [5] and [U a is computed. 
The definition takes into account all experimental un- 
certainties and their correlations, as well as uncertainties 
in the hadronization and underlying event corrections. 
The latter two uncertainties are assumed to be half the 
size of the individual corrections, to be independent of 
each other, and each to be fully correlated over Majct. 
Correlations between the statistical uncertainties are ig- 
nored, since the overlap of the data for the different sce- 
narios is not large. PDF uncertainties are not taken into 
account in the calculations. Otherwise, a theoreti- 
cal prediction affected by large PDF uncertainties and 
in poor agreement with data may get a smaller x^ than 
a prediction with better agreement with data but small 
PDF uncertainties. Therefore, since the PDF uncertain- 
ties are defined differently for different PDF parametriza- 
tions [33, the x^ values would no longer be suited to 
benchmark different PDF parametrizations. This means 
that the x^ values presented here are only a measure 
of the agreement of the central PDF fit results with the 
measured three-jet cross sections. 

The theory results, and therefore the x^ values, depend 
on the choices of as(Mz) and the scales /i^ and /i/ used 
in the computations of the NLO matrix elements and on 
the chosen PDF parametrization. The latter also depend 
implicitly on the value of as{Mz^- AU these dependen- 
cies are shown in Fig. HI where the x^ results are dis- 
played as a function of as(Mz) used in the NLO matrix 
elements and PDFs, for three alternative scale choices 
[ir — — Mo/2 (a), /iQ (b), and 2/xo (c). The results 
for the central as(Mz) choices for the different PDF sets 
are also indicated. For as{Mz) values close to the world 
average of 0.1184± 0.0007 (sij, for all PDF sets, with the 
exception of HERAPDFvl.O, the lowest x^ is obtained 
for the central scale choice /i^ — Pf = /^o- Table [Tl] 
gives a summary of the x^ values, obtained using differ- 
ent PDFs with their default as{Mz) values for the central 
scale fir = M/ = MOj as well as the minimum x^ for all 
scale and as{Mz) choices. From all PDFs, the largest 
X^ values are obtained for CTIO and HERAPDFvl.O 
PDFs, independent of the scale and as{Mz) choices. 
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FIG. 4: (Color online.) The values between theory and data, as a function of the value of as{Mz) used in the matrix 
elements and PDFs. The results are shown for different PDF parametrizations and for different choices of the renormalization 
and factorization scales. The positions of the central as{Mz) values in the different PDF sets are indicated by the markers. 



TABLE II: values between data and theory for different 
PDF parametrizations in the order of decreasing x^, for all 
49 data points. 



PDF set 


Default 
as{Mz) 


X^ at fir = fJ,f = flO 

for default UeiMz) 


Xminimum 


HERAPDFvl.O 


0.1176 


95.1 


81.7 


CTIO 


0.1180 


94.5 


88.2 


ABKM09NLO 


0.1179 


76.5 


76.5 


NNPDFv2.1 


0.1190 


60.5 


59.9 


MSTW2008NLO 


0.1202 


59.5 


59.5 



These are always above > 81.7 for all 49 data points. 
For ABKM09NLO PDFs, which are available only for 
a single value of as{Mz) = 0.1179, the smallest is 
76.5, obtained for = fJ-f = A*o- The best overall 
agreement, corresponding to the lowest values, is ob- 
tained for MSTW2008NLO for the central scale choice 
fir = l^f = A*o and as{Mz) = 0.121 with — 59.5. 
Very close to these are the results for NNPDFv2.1 for 
which the lowest is 59.9 for /Zr = /i/ = /io and 
as{Mz) = 0.123. The large x^ differences between the 
different PDF sets demonstrate the PDF sensitivity of 
the three-jet cross section data. 

In summary, we have presented the first measurement 
of the inclusive three-jet differential cross section as a 
function of M^jct in pp collisions at a center of mass 
energy of ^/s — 1.96 TeV. The three-jet cross section 
is measured in five scenarios, in different rapidity re- 
gions and for different requirements for the jet trans- 
verse momenta. The data are compared to pQCD cal- 
culations in next-to-leading order in the strong coupling 
constant for different PDF parametrizations, by com- 
puting x^ values for different scale choices and different 
as{Mz) values. The best description of the data is ob- 
tained for the MSTW20G8NLO and NNPDFv2.1 PDF 



parametrizations which describe both the normalization 
and the shape of the observed Majct spectra. The PDF 
parametrizations from ABKM09NLO give a reasonable 
description of the data, although with a slightly different 
shape of the Msj^t spectrum. The central results from 
the CTIO and HERAPDFvl.O PDF sets predict a differ- 
ent Majct shape and are in poorer agreement with the 
data. 
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